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With a substantial amount of simulations, we have explored the system across a wide 
range of lattice scales. We have located a lattice artifact, first order bulk transition, have 
studied its properties, and found that the flavor-singlet scalar meson mass vanishes at the 
critical endpoint. We will discuss the lattice phase diagrams and the continuum limits for 
I both a SxSB phase and an IR conformal phase, and compare results with other groups. 
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1. Introduction 



^ It is both pheiiomenologically and theoretically interesting to study strongly coupled 

CN gauge theories with the number of massless fermions, Nf, tuned close to the lower 

^~~^ limit of the conformal window, where the theory remains asymptotically free at 

j^—s high energies while developing a non-trivial infrared (IR) fixed point. While the 

.• system with Nf close to the upper limit of the conformal window can be studied 

(^ perturbatively,^ for systems with Nf close to the lower limit, the behavior is less 

CO clear due to the non-perturbative nature of gauge theories. 

'^ We have been studying the SU(3) gauge theory with 8 and 12 fundamental 

^ fermions on the lattice, using the naive staggered fermion action and the DBW2 

gauge action. We have explored the system through basic low-energy QCD observables 

^ with multiple volumes. In this paper, we will concentrate on our recent findings of 

^ the lattice artifact bulk transition with 12 flavors, discuss the continuum limit, and 

compare results with other groups. We quote results from largest volume available 

at each set of parameters, with finite volume effects of about a few percent. 

2. Lattice critical point 

We have located a line of first order bulk transitional! in the bare lattice parameter 
space: quark mass, niq, and coupling, /3 — 6/g^. This line of first order transitions 
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appears at relatively strong lattice couplings j3 ^ 0.46 and small input quark masses 
TTiq < 0.008. Discontinuities in hadronic observables at the first order transition 
become smaller, as the quark mass, m^, is increased. The system shows a quick 
cross-over region,Elat to, > 0.01, with a marked change in lattice scales, from strong 
couplings to weaker couplings. 
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Fig. 1 . The critical endpoint from extrapolating the mass of the flavor-singlet scalar meson 

The system is approaching a second order transition, or a critical endpoint, as the 
quark mass, mq, is increased while the lattice coupling, /3, is tuned toward the end of 
the first order bulk transition. We have measured the mass of the flavor-singlet scalar 
meson, a, using metastable ensembles at two sets of input parameters, detailed in 
Ref. [2] The mass of the flavor-singlet scalar meson, nia, depicted in Fig. [l] decreases 
on the line of the flrst order bulk transition, as rriq increases, and becomes even 
smaller than pion. To flt the critical behavior, we have used the mean-field scaling 
law, rricr oc {iriq — mr^)^", which happens to describe our data quite well. The fitted 
location of the endpoint from ordered-start configurations (weak coupling) agrees 
with the fitted result from disordered-start configurations (strong coupling). 

This critical behavior of m„ , resembles the behavior found in finite temperature 
transitions.'^ However, our findings with 12 quark flavors at zero temperature occur 
at finite rUq and /?. While the continuum physics dictated by asymptotic freedom, 
can only be extracted at the limit of iriq —> and /? — ?► oo, this critical point certainly 
does not have such a conventional continuum limit, and is irrelevant to the continuum 
physics we are after. Nonetheless, this critical point does possess a continuum limit, 
since the a correlation length diverges at this point, while all other hadronic scales 
decouple from the theory. This second order endpoint thus gives a scalar field theory, 
and is likely non-interactingQ 

The first order bulk transition and critical endpoinir' almost surely depends on 
the lattice action, as a few other groups have found bulk behaviors^^^ with different 
lattice actions and bare parameters. In the parameter space we have simulated, 
however, we did not observe the same characteristics of the bulk transition those 
groups advertised. We acknowledge that certain kind of lattice artifacts are highly 



*We don't have hard proof, but it is hard to construct a nontrivial scalar field theory. 
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likely to exist in different lattice actions under different disguise and with different 
degrees of severity. Our findings of the critical behavior of the flavor-singlet scalar 
meson suggest that caution must be exercised in interpreting results from the scalar 
channel on the lattice, especially when asymptotic scaling is not well established. 

3. Continuum limit 

The continuum limit of lattice gauge theories can be analyzed by the renormalization 
group equation of any physical observable, P{a,g^m), which depends on lattice 
spacing, a, bare lattice couplingrlg, and bare lattice input quark mass, rh — ma. 
Near the continuum limit, P should mildly depend on a, 

where /3{g,rh) and '){g^rh) are lattice version of the renormalization group beta 
functions, and the right hand side of the equation is the usual scaling errorFlln an 
asymptotically free theory, to find the lattice beta functions, one would perturbatively 
compute some observable P, and plug the perturbative form back to Eq. ([I]). See 
Ref. |8]for a massless derivation using the quark potential. One will find the beta 
functions coincide with the continuum ones up to the lowest universal orders |j In an 
asymptotically free theory, this procedure is self-consistent, because g — >■ 0, together 
with a fine-tuned m — > 0, leads to a — > according to the beta function calculated 
from this procedure. 

Away from the physical continuum limit, one can supposedly use a form of the 
observable P that contains lattice higher order operators, say P -^ P + Pi, where 
Pi is irrelevant in the usual continuum limit of physical QCD. With this additional 
term Pi, solving Eq. ^ has the possibility to lead to lattice beta functions, f3i{g,rh) 
and 'yi(g,m), that are completely different from the usual continuum QCD, and 
contain lattice artifacts. One can solve these new lattice beta functions for a — > 0, 
and obtain a new continuum limit at g — >■ gc and m — >■ rh , which is different from 
the usual continuum limit of physical QCD at 5 — > and to — > 0. As long as the 
result of observable P + Pi is consistent with calculations in the limit oi g ^ gc 
and TO — )• m^ , such that Eq. (II]) is self-consistent, we arrive at a lattice artifact 
continuum limit. 

The second order critical endpoint we found with 12 flavors at finite lattice 
coupling /3 and input quark mass niq, presented in the previous section, is such a 
lattice artifact continuum limit. The empirical lattice phase diagram is sketched in 
Fig.[2ja), where the axes are lattice bare parameters. We have two data points, which 
are located on the transition line and extended to smaller masses with a dashed line. 



''Not to be confused with renormalization beta functions, we use /3 and g interchangeably as the 

lattice bare coupling for ease of discussion. 

'^The size of the scaling violation depends on the lattice action and observable. 

■^Minus signs in Eq. Jlk ensure this agreement. 
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(a) Empirical diagram with 12 flavors (b) Infinite volume diagram with IR conformal 

Fig. 2. Lattice phase diagrams 



We note two critical points that correspond to the lattice artifact continuum limit 
and the physical QCD continuum limit, respectively. 

To verify whether the 12-flavor theory is in the conformal window or remains 
in the spontaneously chiral symmetry breaking phase, one needs to know how the 
above picture of continuum limits would change if the underlying continuum theory 
has a non-trivial IR fixed point. In the following discussion, we assume the lattice is 
infinitely large, y ^^ cx). Since the weak coupling limit of the theory is asymptotically 
free, the limit with /3 — > oo and vriq — >■ is a valid continuum limit representing the 
true IR conformal continuum theory. It is conceivable that any non-zero quark mass 
term would break the conformal symmetry and the renormalized theory would lose 
the IR fixed point. If we keep niq = 0, and vary /? away from the weak coupling 
limit, since the renormalized coupling remains the fixed-point value, ^ir, the physics 
governed by the IR fixed point does not change, and neither do the lattice beta 
functions. Thus the theory remains IR conformal with a vanishing lattice spacing, 
a = 0. In terms of Eq. (fTl), any low-energy observable P remains under control of 
the IR fixed point, and is independent of the bare coupling, g, such that increasing 
g from does not need to be compensated by increasing lattice spacing, a. In this 
case, Eq. (jTl) is self-consistently satisfied with a = for m = and < g < 51R. In 
conclusion, shown in Fig.^b), in the infinite volume limit, any point along the line 
with niq = and j3 > 6/5^^ is a lattice critical point, and shares the same continuum 
limit as the weak coupling limit, which describes the same conformal theory at low 
energies. On the other hand, if mq 7^ 0, the theory is no longer conformal, and lattice 
artifacts are present with a 7^ 0, which can be controlled by the bare coupling, /3. 
This proposed scenario appears consistent with the data and we are investigating it 
further. 



4. Mass ratios and universality 

The dimensionless quantity, 7X1^^/ f^^, is shown in Fig. |3] versus the bare input quark 
mass, for both 8 flavors (left panel) and 12 flavors (right panel), respectively. We 
have drawn simple polynomial fit curves on top of our data. The discontinuity at 
0.46 with 12 flavors represents the bulk transition. No bulk 



ruq = 0.006 and /3 
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(a) 8 flavors (b) 12 flavors 

Fig. 3. Scaling of constant physics by vn-^ j f-j^ 

transition has been observed within our simulation parameters with 8 flavors. 

We observed that, at strong couplings, /3 = 0.54 and /3 = 0.45 for 8 and 12 
flavors respectively, the value of m-^/ f^^ decreases as ruq decreases. This behavior 
is expected in a physical QCD simulation, where chiral symmetry is spontaneously 
broken at rUq — > 0, and m,r vanishes while f.^ remains finite in the limit. On the 
contrary, at weaker couplings, /3 > 0.56 for 8 flavors, and (3 > 0.46 after the bulk 
transition for 12 flavors, the m^//7r increases as uig decreases. With 8 flavors, rriTr/ f^ 
is almost constant, and the fitted curves start to bend towards zero at small mq, 
while the fitted curves with 12 flavors continues to grow in the parameter range we 
have simulated. Despite this difference, at weak couplings, TOtt/Ztt with 12 flavors 
resembles that with 8 flavors at relatively large niq > 0.015. 

If we look at the constant physics naively with fixed m-^/fT^ at weak couplings 
of the 12-flavor simulations, in order to reach the massless quark limit, we have to 
decrease /3 and go to stronger couplings. This is in contrast with the continuum 
limit discussed in the last section, and thus we cannot proceed with this procedure. 
The reason for this could be: a) lattice artifacts; b) physical behaviors. 
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Fig. 4. Diagram of dimensionless ratios versus typical lattice correlation length 



We have assembled values of both m^/ f^^ and mT^/mp from our 12-flavor simula- 
tions and plotted them against a typical lattice correlation length l/rrip in Fig. [41 
In the key of the figure, for ensembles with fixed /3 across both sides of the first 
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order bulk transition, weak coupling side (ordered-started) and strong coupling side 
(disordered-started) are labeled with O and D, respectively. The critical endpoint 
is located near the data points with /3 — 0.4626. In the figure, the length of the 
vertical line attached to each symbol represents ?7i^//^, and its position on the 
bottom layer shows the value of niT^/mp versus the corresponding l/rup. Dashed lines 
are polynomial extrapolations inspired by spontaneous-chiral-symmetry-breaking 
(SxSB) behaviors, 



ml = corrigil + cirUg) , 

T-p or fyr =bo+ blTTlq + 63™, 



(2) 
(3) 



acted on each individual channel separately. Assuming S^SB behaviors, l/rUp needs 
to be increasingly larger, with the coupling weakening, to observe the ratios decrease 
as a result of vanishing TOtt. On the other hand, if the theory does indeed have an IR 
fixed point, ratios in the infinite volume limit should approach constant values with 
different couplings, as l/nip — >■ oo, while all the dimensionful quantities vanish in the 
conformal limit. Therefore, smaller quark masses would have stronger differentiating 
power to tell the true behavior of the system. 
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Fig. 5. Comparing dimensionless ratios with results from Fodor et al.^ and LatKMliSI 
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To compare our results with results from other groups, we have included recent 
12-flavor result^ from Fodor et al.^ and LatKMI,'i^in Figure 5 Although results 
of different groups are from different lattice actions, there is the universal behavior 
to be learned from the similar behaviors within the parameter space that has been 
explored. Apparently, both ratios, 111^^/ f.^; and vHt^JvUp are not changing noticeably 
as \/mp increases with fixed /3 within the parameter space limited by the computing 
resources. However, this behavior is not different from the behavior of rrij;/ f-j^, shown 
in Fig. Isla), for a system with 8 flavors. On the other hand, the ratios do seem 
to approach a universal continuum limit at 1/rn.p — >■ cx), if the tendency continues 
without alteration in the region with small quark masses. 



"Results from the largest volume available are taken, and the finite volume effects are small. 
'Values of /,r from LatKMI arc normalized by dividing \/2- 
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5. Summary 

We have located a lattice artifact, critical endpoint, where the scalar correlation 
length diverges, for 12-flavor QCD with naive staggered fermions and the DBW2 
gauge action. In the weak coupling side of the lattice artifact bulk transition, we see 
confined, massive particle states for all simulations done to date. As the input quark 
mass decreases at a fixed lattice coupling, the hadronic scales {rup and /^) decrease 
much more rapidly than the pion mass, such that m^//^ or mT^/mp increases on 
the weak coupling side of the bulk transition. The lattice input quark mass at the 
critical endpoint of the lattice artifact transition gives us a hint of the size of quark 
masses, at that lattice coupling, where the vanishingly small scalar mass can wildly 
distort continuum physics. 

We have argued that an exact same continuum limit of an asymptotically free, 
IR conformal theory can be approached at the massless quark limit and the infinite 
volume limit, with any lattice bare coupling that is weaker than the renormalized 
coupling at the IR fixed point. 

It may not be surprising that much lighter quark masses are required in lattice 
simulations for the SxSB behavior to be numerically indisputable. Nevertheless, how 
to effectively scale towards the continuum limit within our current computational 
resources remains a question for both SxSB and IR conformal scenarios. 
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